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Preface to Third Yearly Report under Contract NAS2-7613 


Work under Contract NAS2-7613 started on July 1, 1973. The 
research goals for a 3 year period stated in tills contract arer 

(a) Assess analytically the effects of fuselage motions on stability 
and random re^onse. The problem is to develop an adequate but 
not overly complex flight dynamics analytical model and to study 
the effects of structural and electronic feedback, particularly 
for hinge less rotors. 

(b) Study by computer and hardware experiments the feasibility of 
adequate perturbation models from non-linear trim conditions. 

The problem is to extract an adequate linear perturbation model 
for the piirpose of stability and random motion studies. The 
extraction is to be performed on the basis of transient responses 
obtained either by computed time histories or by model tests. 

(c) Extend the experimental methods to assess rotor wake-blade 
interactions by using a 4-bladed rotor model with the capability 
of progressing and regressing blade pitch excitation (cyclic 
pitch stirring), by using a 4-bladed rotor model with hub tilt 
stirring, and by testing rotor models in sinusoidal up or side flow. 
Seven reports on the work under Contract NAS2-7613 have been 

submitted, references 1 to 7. Reference 1 completes research goal (a). 
References 2, 4, 6, 7 pertain to research goal (b). It is incomplete 
to date. References 3 and 5 pertain to research goal (c). It also is 
incomplete to date. Reference 3 presents the results of extensive 
frequency response tests* Reference 5 presents dynamic downwash 
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measurements in hovering during harmonic rotor excitation • Three 
manuscripts for publication have been prepared. The first has been 
published as reference 8. Two further manuscripts covering part of 
the material in references 3 and 7 have been submitted to journals. 

The extensive rotor state and parameter identification work with 
coapntev simulated transients was found to be very useful for the 
subsequent data processing of the measured transients. It allowed to 
sort out possible inadequacies of the applied identification algorithm 
and of the inputs from possible inadequacies of the measurements and 
of the applied mathematical rotor model. Rotor state and parameter 
identifications from meastired transients are presently complete for 
hovering conditions using cyclic pitch stirring transients. They are 
presented in reference 7. The corresponding work for forward flight 
conditions and for hub tilt stirring is planned to be completed in 
FY 1977 during an authorized extension of the research contract. It 
is also planned to refine the analytical rotor model utsed for the 
state and parameter identifications to include blade flexibility. 

Rotor state and parameter identification from transients is still 
a field where little experience is available. As elaborated in 
Chapter x of reference 4 it takes four important ingredients to perform 
a successful state and parameter identification from transients; a 
suitable input, a suitable instrumentation for measuring key state 
variables, an adequate mathematical model of the system, and an 
efficient criterion function for the estimation algorithm. In all 
four respects considerable work had to be done in order to finally 
establish a combination of these four ingredients that led to success. 
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The lessons learned in this effort are h||§ed to pave the way for a 
wider application of rotor dynamic perturbation state and parameter 
identifications from transients about non-linear trim conditions, 
both in rotor wind tunnel testing and in rotor fli^t testing. 
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Rotor Dynamic State and Pars?geter Identification From 
Simulated Forward Flight Transients 

K. H. Hohenemser, D. Banerjee, and S. K. Yin 


Abstract 

State and parameter identifications from simulated fozward flight 
blade flapping measurements are presented. The transients are 
excited by progressing cyclic pitch stirring or by hub stirring with 
constant stirring acceleration. Rotor dynamic inflow models of 
varying degree of sophistication are used from a one peurameter inflow 
model (equivalent Lock nunfcer) to an eight parameter inflow model. 

The maximum likelihood method with assumed fixed measurement error 
covariance matrix is applied. The rotor system equations for both 
^fixed hub and tilting hub are given. The identified models are 
verified by comparing true responses with predicted responses. An 
optimum utilization of the simulated measurement data can be defined. 
From the numerical results it can be anticipated that brief periods of 
either accelerated cyclic pitch stirring or of hub stirring are 
sufficient to extract with adequate accuracy up to 8 rotor dynamic 
inflow parameters plus the blade Lock nximber from the transients. 
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Notation 


B 

F. f’ 
Gv» Gu 


J 
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M, Nv, 


Nu 


M 

Pi. ••• P 5 

P 

R 

a 

Cj. 

Cl 

t 

u 

V 


X 

y 

<*k 


blade tip loss factor 

rotor state matrices of size 8x8 and 11 x 11 
respectively 

matrices relating induced flow and coitrol 
variables respectively to the rates of states 

criterion function 

rotor induced flow gain matrix 

matrices relating the state, induced flow and 
control variables respectively to the rotor 
thrust and moment coefficients 

Information matrix 

rotor induced flow parameters 

blade natural flapping frequency 

■Kasurement error covariance matrix 

blade section lift slope 

rotor thrust coefficient, positive up 

rotor pitching moment coefficient, positive nose-up 

rotor rolling moment coefficient, positive to right 

non-dimensional time (period of revolution 2ir) 

time it takes kth blade to move from rear position to 
present position 

control vector 

total mean rotor flow velocity (non-dimensional), or 
noise vector 

rotor state vector 

rotor measurement vector 

hub tilting angle at kth blade, positive up 
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Notation (continued) 
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Superscripts 


T 

it 

Subscripts 

m 

£ 


nose down hxib tilting angle 

left hub tilting angle 

flapping angle of kth blade 9 positive up 

nose down cyclic flapping angle 

left <^clic flapping angle 

blade differential coning angle 

blade Lock number 

pitch angle of kth blade » positive nose up 
nose down cyclic pitch angle 
left cyclic pitch angle 
collective pitch angle 
rotor advance ratio 
total mean rotor inflow velocity 
rotor mean induced velocity 
mean perturbation induced velocity 


non-dimensional 
positive down 


perturbation cyclic induced velocity at blade tip, 
down at rear 

perturbation cyclic induced velocity at blade tip, 
down at right 

rotor solidity ratio, or standard deviation 


time differentiation 
estimated value 
transposed matrix 
equivalent 


measured variable 


empirical 
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Introduction 

This report covers the extension of computer simulation work pre- 
sented in references 2 and 4. Reference 2 refers to rotor d 3 niainic 
state and parameter identifications in forward fll^t (.4 advance ratio) 
using the concept of an equivalent Lock number to appTOxlmate rotor 
dynamic Inflow effects. Two analytical models were tised; a complete 
single blade representation in the rotating frame of reference, and a 
simplified multiblade representation omitting periodic terms and 
omitting multiblade accelerations. Two parameters were identified from 
the simulated noise polluted blade flapping measurements; the equivalent 
Lock number, and the collective pitch emgle. Two transient inputs were 
studied; a rectangular normal flow pulse, and a wave shaped normal flow 
pulse. The simulated measurements with computer generated noise 
pollution were first preprocessed by either a digital filter that took 
out the high frequency noise or by a Kalman filter that used estimates 
of the unknown parameters. The parameters were identified by integrating 
a set of differential equations that sequentially minimized the system 
equation error. The unknown parameters could be particularly well 
retrieved for the wave shaped normal flow pulse. This was true both for 
the single blade model with periodic coefficients and for the approximate 
multiblade model with omitted periodic terms and accelerations. 

Reference 4, Chapter 2, refers to rotor state and parameter 
identifications both in forward flight and in hovering using cyclic 
pitch stirring transients as inputs. The forward flight model was 
limited to the concept of equivalent Lock number, either in the form of 
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a single blade representation vi\ii periodic coefficients, or in the 
form of a multiblade representaticm with constant coefficients. The 
simulated measurements were again polluted by computer generated 
noise. The method used in reference 2 was modified in two ways; first, 
instead of the sequential equation error minimization, a global 
equation error minimization was used that is computationally more 
efficient; second, the Kalman filter to preprocess the measurements 
was used in an iterative way, being ijq>dated whenever new paraiieter 
estimates were available. The iterated equation error estimation with 
updated Kalman filter was compared with a version of the so called 
maximum likelihood methrd where system noise is not modeled . It was 
found that despite somewhat more computer CPU time per iteration, the 
maximum likelihood method was superior because of more rapid convergence 
and because of more meaningful parameter covariances. 

The maximum likelihood method was then applied to th« problem of 
rotor dynamic state and parameter identification from cyclic pitch 
stirring transients in hovering using a time delayed rotor inflow. 

Now 3 parameters were assumed to be unknown; the blade Lock number, the . 
inflow gain, and the inflow time constant. Two inputs were assumed. 

At time zero a step input in cyclic pitch, and at time t = 70 the 
beginning of cyclic pitch stirring with constant stirring acceleration. 
The identification process was started at t = 70. The transient from 
the step input had not completely subsided when cyclic pitch stirring 
began. It was found necessary to include the initial values of the 
blade flapping deflection and of the rotor inflow as unknowns in the 
parameter identification scheme, leading to 7 unknown parameters. 
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Though this required 6.3 CPU seconds per iteration^ (IBM 360/65 computer) 
the second iteration was found to be almost converged, so that the 
total computer effort was moderate. For the tests there is no step 
input of cyclic pitch, so that initial valtie identifications can be 
omitted. 

Reference 4, Chapter 3 develops a method of estimating the 
optimum transient data length that proved to be later very useful in 
all state and parameter identifications. In the maximum likelihood 
method one needs the inverse of the so-called information matrix that 
theoretically gives lower bounds to the parameter covariances. A 
differential equation for the inverted information matrix was derived 
and it was found that its integration results in the definition of a 
data length beyond which no improvement to the parameter covariances 
can be expected. If less data are used in the identification, a rapid 
increase in the parametei^ covariances occurs. Thus a rational way was 
found to avoid errors from insufficient data length and to avoid 
unnecessary computer effort and a degradation of accuracy from 
meaningless additional data. 

As mentioned before, the forward flight studies in references 2 
and 4 were limited to the concept of equivalent Lock number. What 
remains is to apply the concept of a time delayed rotor inflow to 
forward flight conditions and to also consider rotor hub stirring 
transients mentioned in research goal (c). Computer simulations were 
performed to study the feasibility of these types of application of 
the maximum likelihood method. The present report describes this work. 
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Single Blade and Mult*.blade Coordinates 

As before, the nunteical analysis Is performed for an advance 
ratio of .4. Since now dyneunic rotor Inflow that couples the motions 
of the various blades Is Included, a multiblade representation is 
called for. The relations between single blade and multiblade 
variables are: 

Blade flapping angle: = $o + ^II '^k 

Blade pitch angle: ®k * ®o " ®I ‘’^k '^k 

Induced flow: = Vq + Vj(r/R)cos tj^ + vjj(r/R)sln tj^ 

The subscript I refers to forward cyclic flapping, cyclic pitch and 
cyclic induced flow (the inflow is dcwn in the rear), the subscript II 
refers to left cyclic flapping, cyclic pitch and cyclic induced flow 
(the inflow is down to the right). represents differential coning 
for the 4-bladed rotor, whereby one pair of opposing blades cones up, 
the other pair cones down. Though a linear distribution of the induced 
flow over the radius is defined in Eq. (1), this assumption is not 
required for the parameter identification process. Different inflow 
distributions merely produce different values in the identified 
parameters but do not change the form of the equations. 

For the tests the constant collective blade pitch 6o is known. 
Also known are the cyclic blade pitch variables 6j(t)^ ®II^"'^^ 
function of time. The only measured state variables are the 4 blade 
flapping angles B^n,* ®2m» ®3m» ®4m* derived from Eq. (1) 

the multiblade flapping variables are related to the single blade 
flapping variable- by the transformation 
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For the conputer simulations it was assumed that there is additive noise 
in the single blade flapping measurements • Thus the measurement 
equations are 
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The measurement noise covariance R was assumed to be given and not 
changea during the iterations. The innovation is given by (see Chapter 
2 of reference 4) 
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In order to determine the eetimstes 821 03* 64 from the estimates 
for the aultiblade coordinates Bi» Bjj, Bd one needs Eq. (3). 

For the state and parameter identifications from test data a 
subtly different procedure was selected, as discussed in reference 7. 
The measurements 8^, 82^, 8^, 8^^ were first trensfonBed to 
multiblade coordinates by using Eq. (2). These mult ib lade variables 

^om* ^Im* ^Ilm* ^dm considered in the estimation algorithm 

as measured quantities with additive noise: 
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The innovation vector is then given by 


/V 

8 - B 

®dm ~ 


(7) 


The measurement error covariance was not considered given but was 
updated in each iteration* The procedure used for the state and 
parameter identificatiox^ from test data has the advantage of saving 
for each iteration the execution of the transforroation Eq. (3). Also 
in this procedure suitable weights are applied to the test data in the 
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form of the measurement error covariance determined from the preceding 
iteration. 


Inflow Model with One Time Constant 

We adopt here the rotor inflow model of references 9 and 10. 
Eq. (33) of reference 10, written in o\ir notation, reads 
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km 0 0 

0 -kj 0 
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( 8 ) 


Rotor thrust and m<»nent coefficients C,p, C^, are from aerodynamic 
contributions only. Lj. is the enqpirical L-matrix defined in reference 
9. The theoretical values of and kj, using potential flow around 

a solid disk are given in reference 10 as k^ = .849, kj = .113. The 
components of the L-matrix as well as k„, and kj will be identified 
from rotor transient tests. Fr<xn momentum theory one obtains 
according to reference -0 


CL] 


-1 


aa 


2v 0 0 

0 -v/2 0 

0 0 -v/2 


(9) 


yi^ + A(X + v) 

(y2 + x2)1/2 


( 10 ) 


where X and v are the trim values, about which the rotor inflow 
perturbations Vq, v-, vjj are taken. Note that an induced flow 
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trim value is only defined with respect to the axial induced flow. 

For the pitch stirring rotor model we have prior to the start of the 
pitch stirring transient a cyclic pitch trim condition equal to the 
amplitude of the cyclic pitch stirring. This anplitude is considered 
to be a linear perturbation variable and not part of a non-linear trim 
condition. 

For stea^ conditions at advance ratio .4 one obtains frcmn 
reference 9 
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( 11 ) 


In reference 9 cyclic induced inflow was assumed with constant dis- 
tribution over the radius. To adjust for the linear distribution 
assumed in reference 10 and in Eq. CL) the empirical L-matrix given in 
Eq. (11) is changed to 


.5 



0 


0 0 
-2.7 -1.3 

1.3 -4.0 


( 12 ) 


After inversion of this matrix, division by kjjj and kj, insertion of 
p = .4 and rounding off the numbers, one obtains Eq. (8) in the form 
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The rotor moment coefficients and are here assigned with their 
usual sign, positive when the moments from the flow on the rotor are 
nose*-up and to the right respectively. 

Three tj^pes of state and parameter identification from computer 
simulated measurements were performed. In the most general case there 
are 8 unknown inflow parameters in Cq. (13): 
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The 5 parameters to P 5 represent inflow gain constants, while Tq, Tj, 
Tjj are time constants. In this form no assumption is made about the 
relations between the 3 time constants . 

In the second case the theoretical relations Xj = Xjj = 7.5 Tq 
from reference 10 are used. TJiis reduces the number of unknown inflow 
parameters from 8 to 6 ; 
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In the third cass only diagonal taz^ ara ratainad in tha infl.oH gain 
matrix so that P4 P3 * 0. This further raduoaa the number of 
unknown paranetara from 6 to 4. The inflow gain parameters and the 
Inflow time constants will depend on the trim condition and on idte 
rotor advance ratio y . 

Rotor System Equations, Fixed Hub 

The rotor system equatiois are wxdtten in the form 

x = FJi + 6yV + (^u (16) 


with the state vector x given by 
* * *-®o^o®I^I^II^d®d^ 


with the Induced flow vector given by 


Cvoviviil 


and with the control vector given by u'*’ = C6o®I®II^ 


Finally we have a set of equations that gives the aero^mamic 
thrust and moment coefficients in terms of state vector x, inflow 
vector V and control vector u: 


Cp/(ao/2) 

Cj^/(ao/2) 

Cj^/(ao/2) 


Mx + V + N^ u 


(17) 

(18) 
(19) 


( 20 ) 


The matrices F, 6y, Gu» M, Ny, Nu are respectively of size 8x8, 
8x3, 6 X 3, 3 X 8, 3 x 3, 3 x 3. For moderate advance ratio 
neglecting reversed flow effects and assuming 4 straight constant chord 
blades hinged at the rotor center, these 6 matrixes are: 
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Matrix (M^ = B y‘‘/8. Mg * B'® y/3) 
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When inserting Eq*(8) into Eq. (20) the rotor thrust and moment k i 
ficients are eliminated and one (Stains a systwa equation of the ^rm 

X s X + Gu 

tfhere tiie state vector now includes tiie inflow variables 
X*** = C3oSo8lSiSiiPii8d®d'*o'^l'’lI^ 


( 21 ) 


( 22 ) 


where the state matrix F* is of size 11 x 11, and where the control 
matrix G is of size 11 x 3. Note that in the state and parameter 
identifications only 4 out of the 11 state variables a3?e measured. 

In addition to the inflow parameters the rotor system equations 
have as parameters the advance ratio |i that can be considered known, 
the Lock number y and the tip loss factor B. In principle both y 
and B could be identified from the transient test results. We usually 
set B = .97 and then identified y. The system equations also contain 
the blade natural frequency when rotating P that could be easily 
identified from transient test data. In many identifications we 
assumed P as a known quantity. The assumptions made in deriving the 
system equations are less restrictive than they may appear. Only the 
form of the equations is used in the state and parameter identification 
procedure. All parameters are left open and are adapted to the test 
results in such a way that the quadratic differences between predicted 
and actual measurements are minimized. 

Rotor System Equations, Tilting Hub 

As defined in Research Goal (c) transient rotor testing toward 
establishing dynamic rotor wake-blade interactions is to be performed 
both with cyclic pitch stirring and with hub tilt stirring, A second 
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rotor model with the capability of hub tilt stirring is being built. 

The rotor system equations on whidi the state and parameter identi- 
fications will be based are described in the following. 

The second of Eqs. (1) is now replaced by 

“k ” ®I '^k “II "*^k 

where oj and «];j are the forward and left hub tilting angle 
respectively. The blades are again assumed straight and elasticly 
hinged at the rotor center. The blade flapping angles • • • P4 and 
the associated multiblade angles Pj, Pjj are defined not with respect 
to the hub but rather with respect to the space fixed reference rotor 
plane for zero hub tilting angle. Hub tilting aj and ajj has the 
same effect as cyclic pitch application 0j and Pji, except for the 
elastic blade restraining moments that are opposite those from Pj 
and Pii. Thus the control vector is now given by 

u^ = [0Q oj oij] (24) 

instead of by Eq. (19), and elastic restraint terms must be added to 
the Gy - Matrix, since the elastic pitching and rolling moments are 
now proportional to Pj - Oj and Sji “ an respectively instead of 
proportional to Pj and Pjj. 
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G^ - Matrix for Hub Stirring 
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In all other respects Eqs. (16) to (22) Including the expressions for 
the remaining 5 matrices remain unchanged, though the flapping angles 
are now defined differently as compared to the hiib- fixed case. 

The measurement equations (4) are now different since the blade 
flapping angles are measured with respect to the hub. For the conputer 


simulations we have instead of Eq. (4) the measurement equations: 
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The innovation is given, instead of by Eq. (5), by 
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For the state and parameter Identifications from test data a 
somewhat different procedure will be followed, whereby the single 
blade measurements are first transformed into multiblade measurements: 


om 


&dm 


= (1/4) 


1111 
2cost -2sint -2cost 2sint 
2sint 2cost -28int -2cost 
1 1-1 


a L. 


^®1 

"* ^3)111 


(&4 - 


( 27 ) 


The innovation vector is then defined by 

^om “ 

V = 

®dm “ 

The measurement covariance has been assumed given for the computer 
simulations. It will be updated for each iteration for the state and 
parameter identifications from transient test results, as explained before. 

Numerical Results for Cyclic Pitch Stirring Transients 

For the numerical examples presented in the following we assumed as 
given the product of lift slope and blade solidity ratio aa = 2 tt/ 10, 
the advance ratio y = .4, the tip loss factor B = .97 and the blade 
flapping frequency P = 1.20. The Lock number is assumed given as 
Y = 5.0 for some of the identification runs and assumed unknown in 
other runs. The control varidbies 0j, 0jj are assumed given. In 
previous experie'^ce it was found that a moderate noise pollution of the 
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control variables had no substantial eStmct on )the. 
results of the state and parameter identification. For the cases 
presented here 8^ and 6jj were used in the rotor system equations 
without noise pollution* The 4 single blade simulated flapping 
angle measurements were polluted with computer generated Gaussian 
zero mean noise with a standard deviation of = .1, A variety 

of values for the pitch stirring acceleration were studied. It was 
found that the slower accelerations required a longer time span to 
yield the same accuracy of the identified parameters as the faster 
accelerations. Here only the case of a cyclic pitch stirring 
acceleration of m = -.10/u is presented, using a data length of 
12 time units from the start of the transient. Later it was found 
that this data length does not provide optimal data utilization and 
that for better results a data length of 18 time units should have 
been used for m = ~.10/it. 

In preceding studies it was found thar better rotor inflow 
parameter identification can be achieved for progressing cyclic 
pitch stirring as compared to regressive stirring. The reason for 
this experience probably is that for regressive stirring the blade 
natural frequency is resonance excited and that at the resonance 
frequency the dynamic rotor inflow is theoretically zero. All 
presented cases are for progressing cyclic pitch stirring accelerations. 

It was first attempted to solve the problem in the same way as 
describea in Chapter 2 of reierence 4 by including the initial values 
of the flapping angles and of the inflow variables as further unknowns 
to be identi'^'ied. This a^ids aiiother / unknov/ub in the identification 
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routine and led to difficulties. It was then decided to use in the 
identification algorithm perturbations from the initial conditions 
before on-set of cyclic pitch stirring. Thus the initial vilues for 
cdl state variables are zero and need not be identified. 

rig. 1 shows the cyclic pitch stirring input for a progressing 
stirring acceleration of = -.10/ir. The amplitude of cyclic pitch 
stirring is ^ l.S*>. Since the initial value is zero, the maximum 
excursion of 6jj is from 0 to -3°. Fig. 2 shows the flapping 
response of blade number 1 with the computer generated noise pollution 
corresponding to = .10. The polluted values are used in the 

state and parameter identification. Five different analytical models 
will be discussed here. 

(a) Dynamic inflow neglected, Lock number y = 5. 

(b) Quasistatic dynamic inflow, equivalent Lock number y* identified. 

(c) Diagonal L-matrix and inflow time constant identified, y = 5. 

(d) Full L-matrix and inflow time constant identified, y h 5. 

(e) Full L-matrix and inflow time constant and y identified. 

The inflow model given by Eq. (15) is used. In the following table the 

given values (0 or 5), of the parameters are noted. If no value is 

given, the paraneter is identified. 


Model 

Pi 

P2 

P3 

P4 

P 5 

T 

Y 

a 

0 

0 

u 

' 

0 

0 

0 

5,0 

b 

0 

0 

0 

0 

0 



c 



0 

c 



5.0 

d 




1 



5.0 

ft 

[_ 


1 

. 1 

1 

1 J 

1 : 

i 
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The results of the parameter identifications based on the analytical 
models b, c, d, e are given in the following tfl^les: 

Model (b), 4.1 CPU sec/Iteration 



Y* 

Initial Estimate 

575o 

Iteration 1 

3.87 

2 

4.08 

3 

4.09 

M'^^^/Ident. 

.01 

Value 



Model (c), 9.7 CPU sec/Iteration 



‘■i 

?2 

P5 

T 

Initial Estimate 

' .90 

1.08 

.88 

1.08 

Iteration 

1 

.26 

.55 

-.13 

.38 


2 

.36 

.62 

.11 

.53 


3 

.49 

.85 

. 30 

.68 


4 

.54 

.95 

.36 

.74 


5 

.54 

.96 

.37 

.74 

M'^'^^/Iden. 

Value 

.30 

.17 

.19 

.12 


Model (d), 13.3 CPU sec/Iteration 



Pi 

P2 

P3 

P4 

P5 

T 

True Value 

1.00 

1.20 

-.40 

.40 

“.80 

1.20 

Initial Estimate 

.90 

1.08 

-.44 

.36 

.88 

1.08 

Iteration 1 

.74 

.74 

-.71 

.16 

.61 

1.09 

2 

.83 

.90 

-.72 

.29 

.76 

1.19 

3 

.81 

. 86 

-.68 

.25 

.7H 

1.16 

4 

.81 

.87 

-.69 

,26 

.75 

1.16 

5 

.81 

.87 

-.69 

. '"'c 

.75 

1.16 

M"^''^/Iden. Value 

.57 

.63 

.70 

1.9f 

,77 

.51 
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ORIGINAL PAGE IS Po 



22 


Msdel (e), 16 CPU sec/Iteration 



n 

P2 

P3 

P4 

P5 

X 

Y 

True Value 


1.00 

1.20 

-.40 

.40 

.60 

1.20 

5.00 

Initial Estimate 

.90 

1.08 

/ 

-.'44 

.36 

.88 

1.08 

4.00 

Iteration 

1 

1.01 

.63 

■r.89 

.35 

.59 

1.05 

4.85 


2 

.86 

.95 

-.85 

.38 

.81 

1.22 

4.96 


3 

.79 

.81 

-.66 

.23 

.70 

1.10 

4.97 


4 

.82 

.87 

-.71 

.27 

.75 

1.17 

4.99 


5 

.81 

.85 

-.70 

.27 

.74 

1.14 

4.98 


Value 

.60 

.76 

.76 


.86 

.67 

• 

0 

00 


For models (b) and (c) no true values of the unknown parajnevers are 

available, since these models are different from that which generated 

the simulated measurements. For models (d) and the triie 

para-T:eter values are known. One cannot expect that Pi Pp P 3 x for 

Model (c) are the same as for Model (d) or (e). Apparently the much 

lower value for t in Model (c) has to ma^<e up for the omission of 

the off-diagonal terms in the L-matrix, The parameters of Models (d) 

and (e) are aii.:ost the same. Thus identification of y does not affect 

the accuracy of the other parameters, but it does increase the CPU 

time per iteration from 13.3 to 16 seconds. The relative parameter 

- 1/2 

standard aeviations M /Ident, Value are rather large, indicating 
that not sufficient data length has been used. We will see later 
that a substantial reduction in M is possible for greater data 


length • 
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Viarificatlon by Response Conparlsons 

The Identified models have first been verified by comparing the 
responses for the transients used in the identification. The 
comparison in Fig. 3 is between the exact responses - solid lines in 
Fig. 3 and subsequent figures - with the response from the Identified 
model - dot lines in Fig. 3 and subsequent figures-. Fig. 3a is 
for Model (a) without rotor inflow using the cyclic pitch stirring 
acceleration m = -.10/ir. It is seen that substantial errors of 
Model (a) occur in the 3^ response, smaller errors in the 3^ 
response and insignificant errors in the 3^ and 3xi responses. 

These errors are the basis for the inflow parameter identifications. 

Fig. 3b is for Model (b) using the equivalent Lock nuni>er 
concept. The errors in the 3 q ^^<1 @d responses are now insignificant, 
while the 3| and 3|j responses show relatively small errors. Thus 
the equivalent Lock nuni>er concept appears to be quite useful for 
this case. 

Fig. 3(c) is for Model (c) with given Lock niiniber and 4 identified 
inflow parameters; the diagonal L-matrix and one time constant. All 4 
response variables shew negligible errors. Figs. 3(d) and 3(e) are 
for y aels (d) and (e). Though more inflow parameters have been used, 
the response errors are noticeable though still very small. One can 
then conclude that Model (b) identifying only y* and using 4.1 CPU 
seconds per iteration may be adequate in some cases, that Model (c) 
identifying 4 inflow parameters and using 9.7 CPU seconds per iteration 
is excellent, and that the inclusion of the off-diagonal terms in the 
L-matrix is not necessary in this particular case. A sixth model not 
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considered here would probably be useful in state and parameter 
identifications from test transients » namely Model (c) including 
Y as unknown parameter. The Lock nuid>er was found in all cases 
easy to identify accurately* Theoretical values for y not 
always reliable since they involve the blade elastic mode sh^e, 
the lift sic 3 and the exact blade mass distribution. 

The next step in the verification of the identified models is 
to compare responses in transients that have not been used for the 
parameter identification. Figs. 4 and 5 show the true responses to 
a collective and to a longitudinal cyclic unit step input respectively 
as compared to the responses of Models (a), (b) and (e). Again Model 
(b) is a substantial improvement over Model (a) that omits rotor 
dynamic inflow effects. Model (e) gives almost the correct responses. 
Presumably Model (c) i^ith or without identification of y would show 
the same excellent agreement with the assumed ”true'* model. Thou^ 
the inflow parameters to P 5 and t in the various models deviate 
to a certain extent from the true values, the responses predicted 
by Models (c), (d), (e) are in very good agreement with the responses 
of the true model not only for the cyclic pitch stirring transient 
but also for quite different transients. 

Inflow Model with Three Time Constants 

Parcicularly at high advance ratio the time constants for pitch 
and roll could be different, since different participating air masses 
night be anticipated. It, therefore, appeared of interest to study a 
case where 3 instead of 1 time constant are assumed to be lonknown. 

This is a fur*ther extension of Model (e). The complete L-matrix in 
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addition to y was identified. The paran»ters were denned somewhat 
differently from Eq. Q.5X namely by 


1 

o 


1 

o 

o 

0 


1 

0 

6 


*0 


1 

0 

0 

t 


Op 

'>1 

+ w 

0 Tj 

0 


0 

L 22 

^3 



s 

0 Tj 0 


-Cm 

"II 


0 0 

TII 


0 

^32 

1*33 


M 

1 


0 0 Tii 


-Cl 


(29) 

The following table gives the restilts of 3 iterations. The value of 
the criterion function 

N 

J = 11 (Tl„ - R-^ (Yj„ - if) (30) 

1=1 

is given in the last colunaie Contrary to the preceding cases the 
assumed measurement standard deviation for all simulated measured 
variables is a=,05, (R^= 400 I) instead of a = and the 
time used for the identification is 18 instead of 12, Same as 
before the rotor advance ratio is p = .4^ and the cyclic pitch 
stirring acceleration is w = “.l/ir, taken in the progressing sense. 

Also> same as before ^ the time interval for the numerical integrations 
is .1. 


Model (e) with 3 Time Constants, 45 CPU sec/Iteration 



^•11 

Lj2 

^23 

L,2 

^33 

^0 


^II 

y 

J 

True Value 

1.00 

.42 

-.17 

.08 

.17 

1.18 

8.8 

8.8 

3.2 


Initial Est . 

.83 

.33 

-.14 

.06 

.14 

1.35 

7.0 

9.0 

4.5 

8720 

Iteration 1 

1.16 

.29 

-.06 

.08 

.16 

1.15 

3.8 

8.1 

3.1 

951 

2 

.97 

.35 

-.12 

.06 

.18 

1.42 

4.6 

8.5 

3.2 

716 

3 

1.00 

.37 

-.11 

.06 

.17 

1.43 

4.8 

8.5 

3.2 


Error/Ti'ue Value 

0 

.12 

.35 

.25 

0 

.21 

.45 

.03 

0 


/I dent. Value 

.12 

.08 

.71 

.18 

.11 

.13 

.42 

.1» 

.01 
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In comparison to the case for Model (e) with one time constant there is 

better convergence and there are smaller parameter errors. The relative 

“i /2 

predicted parameter error (M /Identified Value) shows the same 
trend as the actual relative error, y having the best and L 23 and Tj 
the worst accuracy. The predicted parameter errors are much smaller 
than for the case of Model (e) presumdily because of the smaller 
simulated measurement error (o = .05 instead of .1) and the longer 
data length (t = IJ instead of 12). The CPU time per Iteration is 
almost 3 times the value for Model (e) with one time constant. In all 
Idle cases presented here the initial estimates are not drastically 
different from the trxie values of the parameters. It was found that 
very much larger errors in the initial estimates can be tolerated 
without affecting the quality of the convergence or of the final 
estimates . 

A somewhat different presentation for verifying the identified 
model is given for this case as compared to the preceding cases. In 
Figs. 6a to 6d the blade flapping response computed with the identified 
model (solid lines) is compared to the simulated measurements (crosses). 
The good fit of the analytical model is evident. Figs. 7a to 7c show 
the dynamic inflow variables Vq, Vj, vjj computed with the identified 
model. Particularly Vjj shows very substantial fluctuations. 

Optimum Data Utilization 

In reference 4 Chapter 3 a differential equation for the inverted 
information matrix M ^ is developed that allows to compute the 
Crairter-Rao lower bound for the parameter covariances vs. the duration 
of the transient losed in the parameter identification process. For the 


RI'TRODUCIBILri'Y OF i 
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example given in reference 4 it was foimd that there enlists an 
optimum data length beyond which no improvement in the parameter 
errors can be expected. It is of Interest to find out whether or 
not such an optimum data length can also be defined for the much 
more complex cases studied in this report. 

The case presented here refers to Model (e) with one time 
constant, the dynamic rotor inflow being represented by Eqs.(13)and 
Q.S). The advance ratio is again y = .4. The simulated measurement 
standard deviation is a = .05. The transient is a progressing cyclic 
pitch stirring excitation with an acceleration o) = The true 

Lock number is y = 5* Fig* 8 shows the predicted parameter standard 
deviations (Cranm^Rao lower bounds) divided by the parameter values 
vs. the duration of the transient used in the identification process. 
The Lock n\imber y shows the lowest relative error, the parameter 
ultimately the highest. It is remarkable that all 7 parameters reach 
their asymptotic relative error at about the same t = 18. Thiis using 
a transient time of t = 12 as was done here for the identification 
of Model (e) with one time constant (but not of Model (e) with three 
time constants), does not give optimal data utilization. On the other 
hand Fig. 8 shows that extending the data length used for the parameter 
identifications much beyond t = 18 would be wasteful of computer 
time and would not result in better accuracies of the parameter 


estimates. 
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Numerical Results for Hub Stirring Tranel«its 

As an alternative to cyclic pitch stirring transients hub stirring 
transients have been considered, and the rotor equations for this 
case have been presented (Eqs. (23|to(26)j). The question is whether 
the accuracy of the parameter Identifications is affected when using 
h\ib stirring instead of cyclic pitch stirring. For blades articulated 
at the rotor center cyclic pitch stirring and hub stirring are 
identical and lead to identical flapping responses. For rotors with 
off-set hinges or for hingeless rotors there are, however, differences 
in the two modes of transient excitation. 

Model (e) with one time constant and the rotor inflow defined by 
Eqs. (L3)and Q.5)has been assumed for the study. The advance ratio is 
y = .4. The stirring acceleration is m = -.l/v in the progressing 
sense. Contrary to the preceding analysis for Model (e) the time step 
has been increased from At = .1 to .2. The duration of the transient 
has also been increased from t = 12 to 24 so that the computer time 
remains the same. 

Model (e). Hub Stirring, At = .2, t = 24, 16 CPU sec/Iteration 



Pi 

P 2 

P 3 

P 4 

P 5 

T 

Y 

True Value 


1.00 

1.20 

-.40 

.40 

.80 

1.20 

5.00 

Initial Estimate 

.50 

.60 

-.20 

.20 

.40 

1.80 

4.00 

Iteration 

1 

1.14 

.57 

-.50 

.50 

.60 

1.11 

4.81 


2 

.94 

1.04 

-.65 

.60 

.82 

1.37 

5.03 


3 

.90 

1.17 

-.59 

.40 

.81 

1.35 

5.06 


4 

.88 

1.15 

-.56 

.33 

.78 

1.32 

5.06 


5 

.88 

1.16 

-.56 

.33 

.78 

1.32 

5.06 

M~^^^/Ident. 

Value 

.25 

.26 

.30 

.75 

.24 

.18 

.02 
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To compare with the equivalent ^cUe pitch stirring casst an analysis 
of this case was made on the same basis* that is using a time step 
At = .2 and a duration t > 24. 

Model (e), Cyclic Pitch Stirring* At s .2* t > 24* 16 CPU sec/Iteration 




^2 



**5 

t 

Y 

Time Value 


1.00 

1.20 

-.40 

.40 

.60 

1.20 

5.00 

Initial Estimate 

.50 

.60 

-.20 

.20 

.40 

1.60 

4.00 

Iteration 

1 

.96 

.91 

-.43 

.44 

.68 

.94 

4.77 


2 

.72 

1.26 

-.46 

.40 

.86 

1.16 

4.96 


3 

.75 

1.26 

i 

• 

to 

.34 

.84 

1.16 

4.97 


4 

.75 

1.26 

CM 

• 

1 

.34 

.84 

1.17 

4.98 


5 

.75 

1.26 

-.42 

.34 

.64 

1.17 

4.98 

M^^'^^/Ident. 

Value 

.27 

.16 

.27 

.27 

.22 

.26 

.04 


Note that in both cases the erx<ors in the initial estimates are much 
larger than in the preceding cases* yet good (»>nvergence is obtained. 

In comparing the hub stirring with the c:yclic pitch stirring case* there 
are only insignificant differences in the parameter errors and in the 
rates of conversion. The predicted parameter errors are also rather 
close except for a larger predicted error for P^ in the case of hub 
stirring. In both cases the predicted and actual errors for the Lock 
number y are quite small. In conparison to the previous analysis for 
Model (e) with one time constant and At = .1* t s 12* both the predicted 
and actual parameter errors are much smaller despite identical computer 
CPU time. 
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Conclusions 

Rotor stato and parameter identifioatlone at .4 rotor advance ratio 
have been performed baaed on simulated blade flipping measurements for 
an analytical rotor nK>dsl tiiat assumes strad^t blades flexibly hinged 
at the rotor center. Re^rsed flox, stall and conpr^sibllity effects 
were omitted, but periodic coefficients in the equations of motion for 
forward flight conditions were retained. All state and parameter 
identlHcations were performed with brief periods of accelerated cyclic 
pitch or hub stirring. The following conclusions can be drawn from the 
study. 

1. The maximum likelihood method with a fixed error covariance matrix 
is well suited for the problem and gives good convergence in all 
cases. 

2. Using a model with an identified equivalent Lock number substcintially 
improves the predicted flapping responses in con^arison to those 
with neglected rotor d}rnamic inflow, but still leaves some errors in 
the responses. 

3. An identified (hmamic rotor inflow model with 4 parameters that 
include one time constant gives almost perfect response predictions. 

4. The addition of 4 more Inflow parameters including 2 further time 
constants and of the Lock number does not affect the rate of 
convergence or the accuracy of the estimates, though it requires 
much more computer CPU time per Iteration. 

Cyclic pitch stirring and hub stirring transients are equally 
suitable for the parameter identifications. 


5 . 
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6. The accupacles with whidi the various parameters can be identified 
are quite different from each other » the Lock nunber having the 
highest accura<^. Despite unavoidable inaccuracies in some of 
the parameters, even relatively large errors in such parameters 
have little effect on the responses. 

7. There is a clearly defined optimal data length that can be 
computed. More data do not bring in^roved accuracy of the estimates, 
while a shorter data length leads to rapidly increasing errors in 
the estimates. 

8. The time step used for the numerical integrations is not a 
critical quantity. Both At = .1 and .2 appear to be satisfactory. 



32 


1* Hoh»n« 08 or» K* H. and Yin, 8. K., "Natiioda Studiaa Toward Sinpllfted 
8otor>Bod]r Dynamto8''« Part I of Flret Yearly Report under Contract 
NA82-7613, dune 1974. 

2. Hohenenaer* K. H. and Yln> S. K., "Conputer Ei^rinents in Pre- 
paratl«i of Syeten Identification fron Transient Rotor Model 
Tests'', Part II of First Yearly Report under Contract NA82-7613, 

June 1974. 

3. Hohenemser, K. H. and Crews, S. T., "Ej^erimsnts with a Four^Bladed 
Cyclic Pitch Stirring Model Rotor", Part III. ef Flxet Yearly 
Report under Contract NAS2-7613, June 1974. 

4. Hohensnser, K. H., Banerjee, D. , and Yin, S. K. , "Meldiods Studies 
on System Identification Frw Transient tetor Tests", Part I of 
Second Yearly Report under Contract NAS2-7613, J\me 1975. 

5. Hohenemser, K. H. and Craws, S. T. , "Additional Eicperlments with a 
Four-Bladed Cyclic Pitch Stirring Model Rotor". Part II of Second 
Yearly Report under Contract NAS2-7613, June 1975. 

6. Hohenemser, K. H. , Banerjee, D. , and Yin, S. K, , "Rotor P 3 rnamic 
State and Parameter Identification From Simulated Forward Fli^t 
Transients", Part I of Third Yearly Report under Contract 
NAS2-7613, June 1976. 

7. Hohenemser, K. H. and Crews, S. T. , "Rotor Dynamic State and 
Parameter Identification FTom Hovering Transients", Part II of 
Third Yearly Report under Contract NAS2-7613, June 1976. 

8. Hohenemser, K. H. and Yin, S. K. , "On the Use of First Order Rotor 
Dynamics in Multiblade Coordinates", 30th Annual National Forum 
of the American Helicopter Society, M^ 1974, Preprint No. 831. 

9. Ormiston, R. A. and Peters, 0. A., "Hingeless Rotor Response with 
Non- IMi form Inflow and Elastic Blade Bending", J. of Aircraft, 

Vol. 9, No. 10, Oct. 1972, pp. 730-736. 

10. Peters, D. A., "Hingeless Rotor Frequenqr Response with Unsteady 
Inflow", Rotorcraft Dynamics, NASA SP 352, February 1974. 

11. Hohenemser, K. H. and Yin, S. K. , "Some Applications of the Method 
of Multiblade Coordinates", J. American Helicopter Soc., Vol. 17, 

No. 3, July 1972. 



33 


Flgura Captions 

Fig. 1 Cyclic Pitch Stirring Inputs with Zero Initial Values for 
Progressing Stirring Acceleration w * -.10 /it. 

Fig. 2 Single Blade Flapping Response with Sitnulated Measurement 
Errors ( s .i. 

Fig. 3 Cyclic Pitch Stirring Respoase Conparisons . 

Fig. 4 Response Conparisons for Collective Pitch ttoit Step Input. 

Fig. 5 Response Cong>arisons for Longitudinal (^clic Pitch Unit Step 
Input. 

Fig. 6 Identified Flapping Responses Using 3 Time Ccmstants. 

Fig. 7 Identified Rotor Dynamic Inflow Using 3 Time Constants. 

Fig. 8 Relative Parameter Standard Deviation vs. Duration of Transient. 
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Fig, 3a, Model (a) 
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Fig. 3b» Model (b) 
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(c) 




Fig. 3d, Model (d) 
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Fig. 3e, Model (e) 
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Fig, 4a, 0Q s u(t). Model (a) 
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Fig. 4b, 0^ s u(t). Model (b) 
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Fig. 4c, 6o u(t). Model (e) 
















Fig. 5b, 0j s u(t), Model (b) 
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Fig. 5c, 0j = u(t), Model (e) 
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Fig. 6c 
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Fig. 64 
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Fig. 7a 
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Fig. 
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Fig. 7c 
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